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1 Introduction 



1.1 The arrow of time 

One of the most self-evident features of our experience is the disparity between time 
evolution of macroscopic systems toward the future and toward the past. So varied and 
pervasive are time-asymmetric phenomena, that it is a remarkable fact and an extraordinary 
example of scientific unification that the vast majority of them can be described in a single 
statement: entropy increases with time. This law — the second law of thermodynamics 
- is a powerful result, but it raises some troubling questions. Boltzmann's entropy is a 
function of microscopic degrees of freedom, yet the laws that govern microscopic systems 
possess symmetries that reverse the direction of time. Where does the asymmetry in time 
come from? 

The natural place to begin is Boltzmann's H-theorem, which purports to derive the 
second law from classical kinetic theory pp. The status of the theorem is less settled than 
often claimed, because it requires the so-called 'molecular chaos' assumption, doubts about 
whose applicability have not been firmly laid to rest. But it is uncontroversial that one way 
or another the second law follows from microscopic physics and an appropriate probability 
distribution over initial microstates. On the other hand, the existence of a low entropy state 
to begin with is more puzzling. Low entropy initial conditions are unnatural in the sense that 
with any straightforward measure they occupy an exponentially small (and unstable) region 
of phase space. If regarded as fluctuations from equilibrium, they are very unlikelyj^] 

So the central problem of the arrow of time consists in finding a justification for the 
so-called past hypothesis [3] — the assumption that the universe had low entropy at early 
times J3 

1.2 Possible solutions 

At the start of the hot big bang era the universe was very smooth, with density pertur- 
bations Sp/p ~ 10~ 5 . This corresponds to a state of very low entropy, because the generic 
fate of classical matter in a decelerating universe is to form an inhomogeneous configuration 
of black holes, i.e. Sp/p ~ 1. It is occasionally suggested that inflation adequately explains 
this fact, since accelerated expansion smooths out inhomogeneities that would otherwise be 
susceptible to gravitational collapse [H [3] . The reason this argument fails [HI El El E] is that 
inflation itself requires extremely special initial conditions to get going. In the simple case of 

As pointed out in [2J, there is no reason at all to suppose that the universe was ever in equilibrium. So 
it isn't quite right to say that low entropy initial conditions are unlikely, hence the euphemism 'unnatural'. 
Nevertheless we would like an explanation for how low entropy initial conditions might have come about. The 
situation is similar to explaining the number of generations of particles. There is no particular reason to say 
that 3 is unlikely, but it is still a mysterious fact that we would like to derive. 

2 A statistical hypothesis is also required, to ensure that the initial low entropy condition is not among the 
pathological few where the entropy gets lower still as time moves forward. This amounts to postulating a flat 
measure (or just about any smooth measure) on the microstates in the initial low entropy macrostate, which 
is already 'natural'. Our aim is to make the past hypothesis just as natural. 
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a single scalar field in a potential V (</>), the equation of state parameter is given by 



For acceleration, we need w < —1/3, which requires the potential to dominate over 'kinetic' 
energy. Since spatial derivatives contribute to the latter, the initiation of inflation requires an 
expanse of spacetime to contain an extremely smooth scalar field. One can see that inflation 
only makes the problem worse in another way by considering a generic state of the sort likely 
to have come about shortly after reheating in an inflationary scenario — most such states 
do not 'un-reheat' when evolved back in time, establishing the specialness of the initial early 
universe state. Similar conclusions apply to the horizon problem inflation purports to solve- 
getting inflation started requires more tuning than simply giving similar temperatures to the 
various causally disconnected patches at decoupling. 

Despite arguments that no dynamical solution of the problem is possible [10], an inter- 
esting scenario was proposed recently by Carroll [11] . In this approach the question of why 
the entropy of the universe was so low in the past becomes meaningless because there is no 
equilibrium state with maximal entropy, and thus any finite entropy is as (un)natural as any 
other. Here we propose an alternative explanation, where the low entropy arises from a phase 
transition in which the strength of gravity increases. A cosmological time asymmetry arises 
from varying spatial asymptotics that leave the local laws of physics time-reversal invariant. 

2 Entropy and gravity 

Gravitational interactions alter the macroscopic appearance of what we consider low and 
high entropy states. Imagine a weakly interacting gas in a box of fixed finite volume. In the 
absence of gravity the gas spreads out. The equilibrium state (i.e. the state with maximal 
entropy) is homogeneous and smooth, and we would be puzzled to find the system in an 
extremely lumpy configuration. 

Now consider the same system, but take account of gravity and imagine that the size of 
the box exceeds the Jeans' length of the gas. The system is unstable to gravitational collapse 
and the result is that most states evolve to highly inhomogeneous configurations of black 
holes — eventually just a single black holej^] With respect to gravitational interactions, the 
equilibrium high entropy states are the lumpy ones, and we would be puzzled to find the 
system in an extremely smooth state. 

The idea we will explore is a transition between these two cases. If, at early times, gravity 
were sufficiently weak (compared to its strength now) , then the natural configurations would 
be homogeneous. Then, as gravity became increasingly strong, at some point the universe 
would find itself in an apparently unnatural state of low entropy (with respect to dynamics 
that include the newly strong gravitational force), despite starting out in a natural high 
entropy state. 

3 We neglect for now quantum and semi-classical effects like black hole decay since these take place on a 
much longer timescale. If the volume of the system is small enough (or if in a cosmological setting the effective 
Hubble constant is large enough), such effects can also be neglected for the reason that the black holes will 
not fully evaporate but will come to equilibrium with their Hawking radiation. 
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Energy hypersurface after gravity turns on 




Figure 1: Phase space of a gas in a box when gravitational interactions are suddenly turned 
on. 

The motion in phase space would look something like that in figure [TJ whose axes represent 
the positions and velocities of the particles in the box. Before gravity turns on, the system 
is confined to a co-dimension one constant E hypersurface, where E is the conserved energy. 
This hypersurface is coarse-grained according to the macroscopic appearance of its contents. 
Shown is one region containing clumpy inhomogeneous states, and another with smooth 
states. Using a 'natural' phase-space measure, the latter occupy significantly more volume 
than the former, and (again using a natural measure) the vast majority of the lumpy states 
consist of configurations which subsequently evolve to homogeneous states. 

Two things happen when gravity turns on (for ease of discussion, we imagine that gravity 
turns on instantaneously). First a negative potential energy is suddenly added to each particle 
pair. This changes the numerical value of the energy from E to E' , since the positions and 
velocities (and hence kinetic energies) are left momentarily unchanged. Second, the foliation 
of the phase space into constant-energy hypersurfaces is altered, because the potential now 
allows energy changes with the momenta of each particle remaining constant. The result is 
that the system leaves the original constant E hypersurface and goes onto the constant E' 
hypersurface it happens to lie on at the time when gravity turns on. 

The E' hypersurface is also coarse-grained into regions according to macroscopic appear- 
ance of the states, but now the regions containing smooth states are small in comparison with 
those containing lumpy states. The system, since it came from a large entropy region of E, 
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necessarily finds itself starting in a tiny smooth region after jumping to E', thus explaining 
the smoothness of initial conditions in a natural fashion. In essence, as gravity is turned on, 
the system is forced from a region of large entropy (with respect to the non-gravitational 
dynamics) into a region with low entropy (with respect to gravity). 

In this simple example, we imagined that the change in the strength of gravity is accom- 
plished by some external agent. If we want the strength of gravity to increase dynamically, 
we need to confront the question of whether the initial state of the dynamical degrees of 
freedom in the gravitational interaction itself must be fine-tuned to allow the transition to 
happen. However, if the time variation in the strength of gravity is truly external, there is 
no such issue; this amounts to solving the arrow of time problem in a rather straightforward 
way - making the laws of physics explicitly time dependent. 

The remainder of this paper develops a simple model that realizes this scenario cos- 
mologically, with one important feature that allows the laws of physics to remain locally 
time-reversal invariant (as deduced from laboratory experiments). The explicit time depen- 
dence will come not from the lagrangian, but rather from spatial boundary conditions which 
we impose on the solutions via a choice of vacuum. The laws of physics are thus locally time 
reversal invariant, but not so on cosmological scales large enough to be sensitive to the effects 
of the spatial boundary conditions. 



3 Review of spatial boundary conditions 

Within the framework of effective field theory, there are three classes of parameters which 
go into generating predictions: 

• the parameters in the lagrangian, 

• the parameters describing the spatial asymptotics obeyed by the fields, i.e. the choice 
of vacuum, 

• the parameters describing the initial conditions of the fields, i.e. the choice of state. 

The lagrangian specifies the equations of motion obeyed by the degrees of freedom, and 
for the field theories that describe our universe well at accessible energies, these equations are 
hyperbolic. Hyperbolic equations generally require two kinds of boundary data to guarantee 
a unique solution. The first is a choice of spatial asymptotics: unless the spatial sections are 
compact, the solution should approach some specified, possibly time-dependent, background 
solution at spatial infinity. This corresponds to a choice of vacuum. The second is a choice of 
field values and momenta over a space-like Cauchy surface, whose spatial values are consistent 
with the choice of vacuum. These are the initial conditions, which correspond to a choice of 
state within the Fock space built over that vacuunj^J 

A field with a given choice of vacuum, i.e. of spatial asymptotics, cannot dynamically 
evolve into a field with a different choicej^] indeed, this is part of the fixing of boundary 

4 Of course, in reality there is no hope of knowing the initial conditions precisely — what one works with 
is a probability distribution over initial conditions. 

5 Solutions with different boundary conditions are separated from one another by infinite energy barriers. 
Quantum mechanical tunneling between vacua is possible via bubble nucleation, but the bubble forms in a 
background of another vacuum whose asymptotics are unchanged. 
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conditions. The space of all field histories is split into superselection sectors according to the 
vacuum to which they belong. About each vacuum, we expand the field in modes. Expanding 
the lagrangian and equations of motion around the vacuum will give the dynamics of these 
mo des0 

The dynamics of the modes will be generated by some hamiltonian. For example, for the 
gravitational modes in GR it is the ADM energy [T2] , 

E = -r- I N ( (2)k " (2) k °) ~ N a Pab r b . (3.1) 

Here the integral is over a 2-sphere at spatial infinity bordering a spatial slice at some fixed 
time, ( 2 )fc is the trace of the extrinsic curvature of the 2-sphere as embedded in the spatial 
slice, and ^ko is the same quantity, though calculated with respect to the background metric. 
N a is the lapse, N the shift, p a b the canonical momentum to the induced metric on the spatial 
slice, and r a the outward-pointing normal vector to the 2-sphere tangent to the space-like 
surface. 

Energy will only be conserved if the background is time translation invariant. In the 
case of a non-trivial, time-dependent background, such as an FRW cosmology, ^ko will have 
explicit time dependence, and so the hamiltonian will depend explicitly on time, and will 
not generate a symmetry. Even though the equations of motion, and hence the local laws of 
physics, are time translation (and time reversal) invariant, the spatial boundary conditions 
break this symmetry. 

We are proposing that spatial boundary conditions represent the 'external force' which 
breaks time reversal invariance in the universe and gives insight into the dynamics that makes 
smooth early conditions more natural. Although a standard FRW universe also breaks time 
reversal invariance in this way, the time-dependence of the hamiltonian is not harnessed 
in a way as to explain the smooth initial conditions; the mere presence of explicit time- 
dependence in the metric does not by itself guarantee that the universe will start out with 
the homogeneous initial conditions that gave rise to our universe 

In what follows, we present a simple model in which there is a background which starts 
in a weak-gravity situation and evolves to a strong-gravity situation. If the spatial boundary 
data are chosen to follow this solution, then we have a situation similar to the gas in the box, 
and we expect that generic initial states will evolve to give a universe like ours. 

4 A model for weakening gravity dynamically 

Our goal is to construct a theory where natural initial conditions lead to a solution in 
which the local strength of gravity evolves dynamically from being very weak in the early 
universe to very strong today. Ideally we would like to have parametric control of the time 
at which the transition happens, and of the ratio of the gravitational couplings at late and 
early times. 

We will accomplish this by building a theory of gravity with two maximally symmetric 
stationary points, one of which has a larger gravitational coupling than the other. The weak 

6 Note that the space of solutions around a given vacuum include all solutions with the given boundary 
data, not just small fluctuations about the vacuum. For example, the asymptotically flat sector of GR should 
include the Schwarzschild black hole. 
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gravity point will be unstable, and initial configurations that start near it will roll down to 
the second (stable) point, in which the local strength of gravity is much stronger. 



4.1 Brans— Dicke theory 

A simple strategy for arriving at a theory of the kind we need is to promote Newton's 
constant to a dynamical variable, and give it a suitable potential. This leads us to consider 
a scalar-tensor action: 

S = Jd 4 x^g (f(cf>)R - l -g{$)g^d^d v <$> - U(<t>)) + C M ■ (4.1) 

Given a vacuum solution in which the scalar takes the VEV 0o> the effective Planck length 
K e s is the coefficient of the Einstein- Hilbert term 

\ = f(<h)- (4.2) 



2K cff 



Two alternative presentations of (4.1) will be helpful for the subsequent analysis. The first 



makes use of the fact that a subset of Brans-Dicke theories can be written as F(R) gravity: 

J d A x V=gF(R) + C m (g»v). (4.3) 

Existence of the two stationary points, and the stability of one of them, as well as the 
requirement that we get different effective Planck constants in the different vacua, then put 
constraints on the function F(R), which we'll find can be met by a simple polynomial whose 
coefficients control the ratio of Planck masses between solutions. 

To see that F(R) actions can be reformulated as Brans-Dicke theories, consider the 
following action^] 

J <Px^-g [F{4>) + F'(<j>)(R - <j>)] + C m (g^). (4.4) 

The equation of motion for the scalar is 

F"{<j>){R-<j>) = Q (4.5) 

which implies R = (ft provided F" ^ oj^] Plugging this back into the action recovers the 
original F{R) action, so the two are classically equivalent^] . A stationary solution in the 
F(R) theory is a constant curvature solution R = Rq to the equation F'(R)R — 2F(R) = 0. 
The 4> equation of motion sets <j> = R, so (f> takes the constant value 4>o = Rq at the stationary 



point. From the Brans-Dicke action (4.4) we then read off the effective value of the Planck 
mass 

' F'( ( f > o) = F'(Ro). (4.6) 



2K eff 



7 See section 5 of [H] for more about scalar-tensor/7 7 '(_R) equivalence and references to the literature, and 
|14| for more on engineering theories with given vacua and Planck masses. 
8 Note that the scalar 4> has mass dimension 2. 

9 F(R) theory is lj — Brans-Dicke theory with a scalar field F'(cj>) and a potential. It is the potential 
that allows some F(R) theories to evade solar system tests. Without it, solar system constraints require 
uj > 40, 000. 
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So far we have presented the theory in two ways: an F(R) action, and an equivalent 
Brans-Dicke theory. It will be useful to consider one further action, found by performing a 
conformal transformation on the Brans-Dicke action to put the Einstein-Hilbert term into 
canonical form. As long as F'((j>) 7^ 0, one can make a conformal transformation to a new 
metric g^ u 

= 2k 2 F\^ v . (4.7) 

The (at this point arbitrary) constant n 2 will have to be positive or negative according to 
whether F'(cp) is positive or negative, and separate conformal transformations may have to 
be made for separate regions in field space. 
The resulting theory reads 



d 4 x 



(4.8) 



where 



1 3 F" 2 rhF' — F 

Whether the scalar is normal or ghost-like depends only on the sign of k 2 , which is the same 
as the sign of F' in the region under consideration. 

Once the conformal transformation is made the scalar 4> is minimally coupled to gravity, 
so vacuum solutions and their stability can be straightforwardly read off from the potential 
y (0). One can check that the points <po at which V'(4>o) = correspond to curvatures i?o = </>o 
satisfying F'(Rq)Rq — 2F(Rq) = in the original F(R) theory. The second derivative V"(4>) 
tells us about the stability of our vacua. Evaluated on stationary solutions, we have 

{M = (2*2)2^)3 • (410) 

However, one must keep in mind that in this frame matter is not minimally coupled. Matter 



particles do not follow geodesies of g^ v . It is more natural to think of the Jordan frame (4.1 ) 
as physical. 



To recap, the equivalent forms of the action are (4.3), the F(R) action, (4.4), the scalar- 
tensor form where matter is minimally coupled, and (4.8), the scalar-tensor form with a 
canonical Einstein-Hilbert term. In the first, the simplicity of the F(R) action facilitates the 
search for the members of this class of theories with the properties we want, the second is 
useful for studying the behavior of matter at a fixed <j) stationary point, the third makes the 
dynamics of 4> more transparent. 

Constraints on F(R) 

With the above in mind, we seek a function F(R) that satisfies the following criteria: 

1. F'(R)R - 2F(R) = should have solutions at R = and at R = A, where A > 0, and 
there should be no other solutions between these two. We choose A > so as to allow 
the flat FRW ansatz. We want to start off in the R = A solution and roll down to 
the R = solution. There will be no exponential expansion at the endpoint, and we'll 
expect solutions to approach a standard power law FRW universe. 
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2. We want F'(R) > everywhere between the two solutions. This will ensure that the 
conformal transformation is valid in this entire region. 

3. We want F"(R) > everywhere between the two solutions. This will ensure that the 
equivalence between (4.3) and (4.1) is valid everywhere in the region, and also that the 
effective Planck mass at the high curvature point is larger than at the zero curvature 
point. 

4. We want 

V"(0) > 0, 

so that the zero curvature solution is stable. As a consequence of the other conditions, 
the high curvature solution will automatically be unstable, and since there are no 
stationary points in between, we can expect that any cosmology starting near the high 
curvature solution will dynamically evolve towards the low curvature solution. 

If we try to satisfy these constraints using a polynomial, we find they cannot be satisfied by 
a quadratic, or by any polynomial of odd order. For a quartic, the constraints amount to: 

1 



F(R) = R + CR 2 + 



A 2 



with either: C > and < E < (2A 3 )" 1 , or C = and E 
allow C to vary: 

R 4 



2E\ ! R 3 + ER A 

(2A 3 )- 1 . 



(4.11) 
We will fix E but 



F(R) = R + CR 2 + 



2A 3 ' 



The factor by which G then increases as the universe moves from R = A to R = is: 



G eff {R = 0) = K? eff (R = 0) = F'(X) 
G eff (R = X) K\ ff {R = \)) F'(0) 



3 + 2CA. 



(4.12) 



(4.13) 



4.2 Our model 



With our choice of F(R) (4.12) the action (4.3) becomes 



-9 



1 

2^2 



R + CR 2 + 



R*_ 
2A 3 



The equivalent scalar tensor model is 

d^x^g- [f(<t>)R - U (<(>)] + C m (9^) 



where 



F'( 



1 

2^2 



1 + 2C</> + 



A 3 



^F\4>) - F(0) 



1 

2^2 



C(f> 2 + 



30 4 
2A7 



(4.14) 
(4.15) 

(4.16) 
(4.17) 



The potential in Einstein frame (see figure [2]), from which we can read off the dynamics, is 

4>F' — F A 3 (34> A + 2CA 3 (/> 2 ) 



( 2k 2)2 F /2 4 /c 2((2C0+1)A 3 + 20 3 ) 2 ' 



(4.18) 
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v(4>) 




Figure 2: Potential in Einstein frame. 



Dynamics 



Figures [3}{5] show the time evolution in the model A = 5, C = 2, with initial conditions 
= 5, and (j> given an initial velocity ~ — 1. The Hubble parameter starts out approximately 
constant, as would be expected in an almost stable point with R > 0, and then descends 
toward zero as the spacetime approaches a power law FRW universe. Figure [4] shows the 
scalar rolling towards the origin during the transition, and figure [5] shows the consequent 
behavior of the effective gravitational coupling (normalized so that G pre sent day = 1) which 
increases by a factor of 3 + 2CX = 23. 

The matter content has been set to zero, but the picture remains qualitatively unchanged 
after adding matter. In fact, adding matter prolongs the time spent near the unstable sta- 
tionary point. 

We can gain parametric control over when the transition happens by adding a scalar 
kinetic term to the Jordan frame action, 
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f^R-lh^f-U^) 



(4.19) 



where h is a constant. Adding this term does not alter the location, stability, or physical 
properties of the vacua, but by making h large, we damp the motion of <j), prolonging the 
time it spends near the unstable stationary point 
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The numerical simulations above show that large increases in the value of Newton's con- 



stant can be brought about in the model (4.14). We therefore suggest that a scenario like this 



can naturally account for the observed smoothness of the CMB and the low entropy initial 
conditions of our universe. 



The additional kinetic term also suppresses spatial fluctuations of c 
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Figure 3: Hubble parameter. 
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Figure 4: Scalar VEV. 

There is of course fine tuning in the lagrangian and in the spatial boundary conditions. 
Indeed, F(R) lagrangians are very unnatural from a quantum field theory viewpoint, and 
there is no a priori preferred reason to choose spatial boundary conditions as we have. We 
therefore don't claim to have found a completely natural and tuning-free solution to the 
arrow of time problem, or a complete theory of the early universe, but only the more modest 
goal of giving a concrete model in which the initial conditions can appear highly tuned, but 




11 



G 



1.0 - 

D.8 - 

0.6 - 

0.4 - 

0.2 - 



100 200 300 400 

Figure 5: Effective G. 



aren't. 

5 Discussion and possible issues 

The road to a natural explanation of the low entropy initial conditions of the universe 
is fraught with difficulties, and one must be particularly sensitive to implicit assumptions 
tantamount to the conclusion itself, which thereby render the putative solution question- 
begging. The subtlety of such considerations makes it difficult to confidently declare a given 
scenario to be free of fine-tunings. Several apparently questionable assumptions were made 
in the model above, some of which we discuss here, mindful of the fact that others may be 
lurking undetected. 

The simulations were carried out with a homogeneous scalar field. This serves only to 
identify the background solution setting the spatial asymptotics. One must imagine a set of 
initial fluctuations on top of these solutions. Indeed, if we are allowed to choose homogeneous 
initial conditions then the very problem we are addressing — smooth initial condition - 
evaporates. One might worry that the increase in the strength of gravity disappears with 
more general initial conditions. However, we argue that the homogeneity of 4> is not necessary 
in the above scenario (except for computational convenience). The phenomenon we rely on 
is just the rolling of the expectation value (0), and unlike inflation, this does not require a 
smooth field. Our background solution, which is indeed fine tuned (it starts near the unstable 
point, whereas it could have started anywhere), only serves to set the spatial asymptotics, and 
allows us to argue that generic configurations of initial fluctuations will lead to a universe like 
ours. We do not address how inflation may tie in with the weakening of gravity, or whether 
it can be realized in a model like this, since this would obfuscate our main point. 

Our scenario does pass one crucial test, avoiding the "double standard" [E]. Often 
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natural-seeming scenarios, such as a ball rolling down a hill and coming to a stop due to 
friction, become completely unnatural viewed in time reverse. In our scenario, the time 
reverse is completely natural and goes as follows: gravity is strong and the universe starts 
in a clumpy state. Then at some point gravity turns off. The clumps start to dissolve into 
a smooth state. The state at all times is natural with respect to the laws operating at those 
times. Passing the double standard test requires that the universe spend a sufficient period 
of time in the weak gravity phase, otherwise there will not be sufficient time (in the time 
reversed sense) for a generic initial clumpy state to smooth out. In our model, we accomplish 
this by tuning the parameter h in the lagrangian. As we've mentioned, the parameters in the 
lagrangian are fine-tuned; it is only the initial conditions that are allowed to be generic. 

Again, our model itself is meant only to provide a concrete realization of the underlying 
idea of weakening gravity, not to be a complete theory of the early universe. As such, the 
idea of an early phase characterized by weaker gravity should be kept distinct from the 
particular features of the model presented. Indeed, our particular model (making use of 
spatial boundary conditions as a time dependent driving force) may not be the only or best 
approach to weaken gravity in the early universe. Another class of models that may do the 
job is presented in [16] . The essential, model independent, point is that a gravitational theory 
whose strength was sufficiently weak at early times can make smooth initial data the norm, 
and so potentially offers a different perspective on the the usual puzzle of low entropy initial 
conditions. 
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